NASA Technical Memorandum 101721 


Effects of Simplifying Assumptions 
on Optimal Trajectory Estimation 
for a High-Performance Aircraft 

Lura E. Kern, Steve D. Belle, and Eugene L. Duke 


April 1990 


( NAS A- TH- 101 7 2 1 ) EFFECTS OF SIMPLIFYING 
ASSUMPTIONS ON OPTIMAL TRAJECTORY ESTIMATION 
FOR A HIGH-PERFORMANCE AIRCRAFT (NASA) 

-7 p 


CSCL OIC 


G 3 / 0 8 


N90-25142 


Unci as 
G289157 


NASA 

National Aeronautics and 
Space Administration 


i in wnmiiniiiniiiim m mi m n i mi ni' hi irinivi i f i i s i iiiiiintinvnm ^ m i 


NASA Technical Memorandum 1 01721 


Effects of Simplifying Assumptions 
on Optimal Trajectory Estimation 
for a High-Performance Aircraft 

Lura E. Kern and Eugene L. Duke 

Ames Research Center, Dryden Flight Research Facility, Edwards, California 
Steve D. Belle 

PRC Systems Services, Edwards, California 


IWVSA 

National Aeronautics and 
Space Administration 
Ames Research Center 
Dryden Flight Research Facility 
Edwards, California 93523-0273 



il III llilllillll 111 h Hill III! ill ill I II II llllll II ill II II II II I II II 1 1 II i II II lllliiill III 1 1 



CONTENTS 

INTRODUCTION 1 

NOMENCLATURE 1 

9 

Subscripts 

ENERGY STATE APPROXIMATION 2 

Background 2 

Flight Envelope Calculations 4 

Trajectory Calculations ‘ 

Aircraft Model ® 

MODELING ASSUMPTIONS 6 

Atmospheric Models ® 

Gravitational Model Equations ' 

Equations of Motion ' 

Constraint Conditions ° 

RESULTS 9 

Effects of Atmospheric Models 9 

Effects of Gravitational Models 9 

Effects of Equations of Motion 9 

Effects of Constraint Conditions on Trajectory Calculations 10 

CONCLUSIONS 10 

REFERENCES 11 

FIGURES 11 


PRECEDING PAGE BLANK NOT FILMED 


I II I li dill iiiliillll iliilUllllill 



INTRODUCTION 


When analyzing the performance or an aircraft, certain simplifying assumptions, which decrease the complexity 
of the problem, can often be made. The degree of accuracy required in the solution may determine ihc extent 
to which these simplifying assumptions are incorporated. A complex model may yield more accurate results if it 
describes the real situation more thoroughly. However, a complex model usually involves more computation time, 
makes the analysis more difficult, and often requires more information to do the analysis. Therefore, to choose 
the simplifying assumptions intelligently, it is important to know what effects the assumptions may have on the 
calculated performance of a vehicle. This paper examines several simplifying assumptions, compares the effects of 
simplified models to those of the more complex ones, and draws conclusions about the impact of these assumptions 
on flight envelope generation and optimal trajectory calculation. Models which affect an aircraft are analyzed, but the 
implications of simplifying the model of the aircraft itself arc not studied. The examples in this paper are atmospheric 
models, gravitational models, different models for equations of motion, and constraint conditions. 

The results are calculated using the energy state approximation (Bryson, Dcsai, and Hoffman, 1969). In the 
energy state approximation, the aircraft is modeled as a point mass in the vertical plane, and specific energy is used 
as a state variable. The results arc for a high-performance aircraft in minimum time-to-energy trajectories. 

NOMENCLATURE 


a 

speed of sound, ^ 

D 

drag, lb 

E, 

specific energy, ft 

f 

function, or forces, lb 

G 

ft 3 

gravitational constant, 

9 

acceleration due to gravity, 

h 

geopotential altitude, ft 

J 

cost function 

L 

lift, lb 

Lm 

temperature lapse rate, ^ 

M 

Mach number 

Me 

mass of the Earth, lb 

Mo 

mean molecular weight 

m 

mass of the vehicle, slug 

P 

1 h 

pressure, jp- 

PLA 

power lever angle 

Ps 

specific power, ^ 

9 

dynamic pressure, ^ 

R 

range, ft 

Re 

radius of the Earth, ft 


sec" 



ft IK 

R* universal gas constant, (ib-mol) °R 

T thrust, lb, or temperature, °R 

l time, sec 

v velocity, 

w fuel weight, lb 

w fuel flow, jjj 

z geometric altitude, ft 

a angle of attack, rad 

7 flightpath angle, rad, or ratio of specific heats 

p density of air, 

first derivative 

Subscripts 

ave 
b 

idl 
h 

Urn 
max 
min 
re f 
v 
x 
z 

ENERGY STATE APPROXIMATION 
Background 

The trajectories discussed in this paper are calculated using the energy state approximation theory (Bryson, 
Desai, and Hoffman, 1969). This theory models the aircraft as a point-mass and assumes a flat, nonrotating Earth, 
which is reflected in the equations of motion used. The general equations of motion for a point-mass in the ver- 
tical plane arc derived by Micle (1962). They can be given by, first, summing the forces in die direction of the 
velocity vector 

T x cos a — T 2 sin a — D — mg sin 7 = mV (1) 

where T is thrust, a is angle o f att ack, D is drag, mg is the vehicle’s weight, 7 is the flightpath angle, m is the 
vehicle s mass, and V is the first derivative of velocity. By summing the forces perpendicular to the velocity vector 

T x sin a + T z cos a + L - mg cos 7 = mV ^ (2) 


average 

base 

idle 

in the horizontal plane 
limit 

maximum 

minimum 

reference 

in the vertical plane 
j-body axis of the vehicle 
z-body axis of the vehicle 
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where L is lift, V is velocity, and 7 is the first derivative of the llightpath angle. These equations assume that 
the vehicle thrust vector is fixed with respect to the vehicle body axes and acts in both the i and z body axes of 
the vehicle. 

In the energy state approximation, the total energy per unit mass, or specific energy, ( E „ ) is considered a state 
variable and is derived from the energy equation (Bryson, Desai, and Hoffman, 1969) 


1 2 

E = mgz + -m” 


(3) 


where E is energy, and z is geometric altitude. This is divided by the vehicle weight to give the specific energy 


E a - z + — 

2 ? 

When this quantity is differentiated with respect to time, the equation for specific power (P a ) is given 

„ . vv 

P a = z+ 


(4) 


(5) 


An expression for V can be derived from equation 1 resulting in 


T x cos a — T z sin a - D - mg sin 7 
V ~ m 

(6) 



An expression for z can be derived from the definition of 7 


• -1 * 
7 = sm - 

(7) 

so that 


i = V sin 7 

(8) 

Substituting equations 6 and 8 into equation 5 results in: 


Tr (T% cos a - T z sin a - D - mg sin 7 N 

P, = y Sl n 1 + V(, m 9 ) 

(9) 


By subtracting out the V sin 7 terms this equation gives an expression for P a . 

Tx cos a — Ti sin a — D 

mg 

Equation 10 gives P s as a function of velocity, thrust, angle of attack, and drag. Drag, however, is a function of 
Mach number, altitude, and angle of attack. Thrust is a function of power lever angle (PLA). Therefore, it can be 
seen that P 3 is a function of altitude, Mach number, angle of attack, and PLA and is independent of flightpath angle. 

Given the expression for P s , it is possible to derive the cost function for a minimum time-to-energy path. The 
minimum time-to-energy cost function is defined as (Erzbcrger, Barman, and McLean, 1975) 

J=[ i dt HD 

Jo 

where J is a cost function and t is time. 
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From equation 5, specific power is defined as 


Ps = 


dEs 

dt 


Substituting equation 12 into equation 11 gives 


J = 


f E • dE, 
o P, 


( 12 ) 


(13) 


To find the minimum for equation 13, it is sufficient to minimize the quantity inside the integral ■ Since 
the energy state approximation assumes that cost functions arc computed along curves of constant specific energy, 
it is sufficient to minimize the quantity p-. This is the same as maximizing P„. The cost function for a minimum 
timc-to-cncrgy trajectory, then, can be given by: 


« 


maximize J — P, 


(14) 


Flight Envelope Calculations 

The flight envelope can be calculated by solving the general equ ations of motion (eqs. 1 and 2) for trim in straight 
and level flight. For this ease, 7 = 0, V = 0, and 7 = 0 . Therefore, the general equations of motion for a point 
mass in the vertical plane are reduced to: 


T x cos a — Tz sin a - D = 0 
T x sin a + T z cos a + L - mg = 0 


(15) 

(16) 


For actual implementation in a numerical nonlinear equation solver, the above equations can be written as: 


fh(ot, PLA, 2, M) 


T x cos or — T z sin a 
~~D 


(17) 


fv(ot, PLA, z, M) 


L + T x sin a - T z cos a 
mg 


(18) 


The trim equations of motion, then, are defined by selling f v (oi, PLA, z, M) - 0 and fh(a, PLA, z, M) = 0. Note 
that these equations make no attempt to force the pitching moment equation to zero. While th is is consistent with 
the assumption of the vehicle as a point mass, the effects of the vehicle’s longitudinal surfaces are neglected. 


The flight envelope is calculated in four parts (fig. 1): the angle-of-attack boundary, the thrust boundary, the 
Mach boundary, and the dynamic pressure boundary. The angle-of-attack boundary is defined by setting the angle of 
attack equal to the maximum angle of attack of the vehicle. An altitude is picked, and the trim equations of motion 
arc solved for that altitude. This gives 

fh( “Ktoi PLA, Z re f, M) = 0 (19) 

and 


U(on im ,VLK,z Te} , M) = 0 


( 20 ) 


This gives two equations and two unknown variables, which can be solved using a numerical nonlinear 
equation solver. 


I 
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The thrust boundary is calculated using the same method. The PLA is set to its limiting value and a Mach number 
is chosen for the calculations. The two trim equations of motion arc then solved for the corresponding angle of attack 


and altitude: 


and 


PLAh'mj Z } A f r tf) 0 

PLA/| m , Z y Mref) — 0 


( 21 ) 

( 22 ) 


The Mach boundary is given by setting the Mach number equal to the vehicle’s Mach limit. 

The dynamic pressure boundary is calculated by solving the equation relating dynamic pressure and velocity, 
using the dynamic pressure limit of the vehicle 


V,./ = (23) 

where q is the dynamic pressure and p is the density of air. The Mach number is found from the resulting velocity, 
and the trim equations of motion arc solved using a reference altitude 


and 


A(a,P LA ,Zre/, M re f) = 0 

(24) 

f v (a, PLA, Z Te f, Mref ) = 0 

(25) 


Trajectory Calculations 

To calculate trajectories, the cost function is maximized along curves of constant specific energy within the flight 
envelope of the vehicle. For the numerical results in this paper, this was done in the following way: 

1 . An array of specific energy curves encompassing the vehicle’s flight envelope was determined (fig. 1). 

2. A point along the first curve, representing the lowest value of specific energy, was chosen, yielding an altitude 
and a corresponding Mach number. 

3. This altitude and Mach number point was entered into an optimization routine that varied angle of attack and 
PLA to determine a maximum value for the cost function, specific power, at that point. 

4. A different altitude and Mach number point on the first curve was then chosen, and the optimization routine 
calculated a maximum value for the cost function at the new point. 

5. These two maximum values for the cost function were compared and the higher of the two was stored. 

6. The maximum cost function value was calculated at other points along the first energy curve and the highest 
value found along the curve was stored. This value identified the optimum point on that energy curve. 

7. The optimum point was then calculated along the remaining energy curves within the flight envelope. The 
altitude and Mach number path connecting these points defined the optimal trajectory. An optimal trajectory 
is shown in figure 1. 

Parameters such as the time, range, and fuel used along the trajectory can be calculated between energy curves. 
This is not the same as the analysis at each point along an energy curve, which assumes steady, horizontal flight. In- 
formation for the total trajectory can be gained by summing the time, range, and fuel used between each 
energy curve. 
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The time it takes to go from one energy curve to another is given by 


At = 



( 26 ) 


Range ( R ) is then calculated based on the lime 


A R = V cos 7 A t ( 27) 

where 7 is given in equation 7. Finally, fuel used is a product of the fuel flow (ti;) multiplied by the time increment 

Aw = wAt (28) 


where w is fuel weight. 

Aircraft Model 


As stated previously, the energy slate approximation models the aircraft as a point mass in the vertical plane. 
This paper does not assess the assumptions made about the aircraft, but concentrates instead on assumptions made 
about outside parameters which will affect the trajectory calculations. These tr ajec t ory c alculations were based on a 
mathematical model of a typical modem, supersonic fighter aircraft. This model supplied lift, drag, thrust, and fuel 
flow as a function of Mach number, altitude, angle of attack, and PLA throughout the flight envelope. 

MODELING ASSUMPTIONS 


Atmospheric Models 



Models of the atmosphere arc used to determine altitude-dependent parameters such as air density, pressure, 
temperature, and speed of sound. This paper looks at the 1962 standard atmosphere (NASA/U.S. Air Force/U.S. 
Weather Bureau, 1962), as well as models for a standard day, cold day, and hot day for Edwards Air Force Base, 
California (Johnson, 1975). The equations used to determine atmospheric-dependent parameters are given here. 

The general equation for temperature is (NASA/U.S. Air Force/U.S. Weather Bureau, 1962) 

T=T b + L m (h-h b ) (29) 


Here, L m is the gradient of the temperature with gcopolcntial altitude, h b is the geopotential altitude at the base of 
a particular layer by a spcci fic value of L m , and T b is the value of T at altitude h b . Temperature profiles for the four 
day types arc given in figure 2 . 


The general equation for pressure (P) is dependent on the value of the gradient (L m ). If L m 7 ^ 0, then P is 


given by 


ijfo 

p_ = / n 

P b \T b +L m Hj 


and if L m = 0 , P is given by 



Here, H = h — h b . Mo is the mean molecular weight, and R* is the universal gas constant. 


(30) 

(31) 
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(32) 


Density is calculated from the pressure and temperature: 

Mi 0 P 


P = 


R*T 


and, finally, the speed of sound is calculated from the temperature: 


a i Mo 


(33) 


Gravitational Model Equations 


Acceleration due to gravity can be calculated several ways. The two ways that make sense for a model that 
assumes a flat, nonrotating Earth arc: 1) calculating gravity as a function of distance from the surface of the Earth, 
and 2) giving gravity as a constant value, assuming no variation with altitude. 

The acceleration due to gravity as a function of distance from the Earth is (NASA/U-S. Air Force/U.S. Weather 
Bureau, 1962) 

r' \4 

(34) 


GM t 


9 ' (Re + z) 2 

where G is a gravitational constant, M t is the mass of the Earth, and R c is the radius of the Earth. 

A simpler means of determining the acceleration due to gravity is to assume that it does not vary with distance 
from the Earth. The value for g, then, is given by the constant 


g = 32.174 


ft 

sec 2 


(35) 


Equations of Motion 

The general equations of motion for a point mass in the vertical plane are given by equations 1 and 2. Several 
assumptions can be made which affect the flight envelope and trajectory calculations in two distinct ways. First, 
assumptions for the equations of motion will affect the calculations of the flight envelope. Second, the equations of 
motion are used in calculating specific power, which is used as the cost function for calculating minimum time-to- 
cnergy trajectories. 

To calculate the flight envelope, the aircraft is assumed to have straight-and -level flight at trim conditions, giving 
equations 15 and 16: 

T x cos a — T z sin a — D = 0 
T x sin a + T t cos a + L — mg = 0 

These give the most general method of modeling the trim equations of motion used in calculating the vehicle's 
flight envelope. 

If the thrust vector is assumed to be aligned with the i-body axis of the vehicle, then there is no thrust along the 
z-body axis, and the equations simplify to 

(36) 

(37) 


T cos a - D = 0 
T sin a + L — mg = 0 
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To further simplify the equations, it can be assumed that the angle of attack is small, in which ease cos a ps 1 
and sin a & 0. This gives 

T-D= 0 (38) 


L — mg = 0 


(39) 


These same assumptions can be used in calculating specific power. The most general equation for speci fie power 
is given by equation 10: 

T, cos a — T z sin a — D\ 


Ps=V 


mg 


If it is assumed that the thrust lies along the body axis of the vehicle, then T z is zero and the equation for specific 
power becomes: 

/ r I ' mm I 1 \ 

(40) 


P, = V 


/T cos a — D\ 
\ mg J 


Finally, assuming that the angle of attack is small gives: 

r.-vM 

\ mg J 


(41) 


Constraint Conditions 

Constraint conditions limit the possible solutions to an optimization problem. The more constraints put on the 
problem, the smaller the range of solutions becomes. Because constraint conditions arc applied to the optimization, 
the choice of the constraint conditions affects only the optimal trajectory and has no affect on the flight envelope. 

The simplest ease in this paper is an unconstrained trajectory. This means that there are no constraints on the 
cost function, and only necessary limits arc put on angle of attack and PLA to keep them within realistic values. The 
maximization problem for optimizing unconstrained trajectory is 

maximize /(a, PLA, h, M) 

subject to a m i„ < a < a max 

PLA idl < PLA < PLA™, 

The quantity J(a, PLA, h, M) is the cost function being optimized, specific power. 

Because the energy state approximation assumes that the aircraft is in straight and level flight at each point, this 
can be reflected in the constraint conditions. To do this, a constraint equation that balances the equations of motion 
in the vertical direction is added. The maximization problem for optimizing trajectory constrained to level flight is: 

maximize /(a, PLA, h, M) 

subject to / t (a:, PLA, h, M) = 0 

O^min ^ OC Ol max 

PLA idi < PLA < PLA moi 

Finally, the trajectory can be constrained to realistic flightpalh angles. This can also be thought of as constraining 
the velocity in the vertical direction, so that the rate of climb (z) can never exceed the total velocity. This can be 
seen from the definition of the flightpath angle 

• -i 2 
T, = sm - 
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The maximization problem for optimizing a trajectory constrained to realistic flightpath angles is 

maximize /(or, PLA,/i, M) 
subject to /„( a, PLA, h, M) = 0 
-7 <T<I 

C^rntn ^ ® ®mai 

PLA, d; < PLA < PLA rooi 


RESULTS 

Effects of Atmospheric Models 

The effects of the atmospheric models on the trajectory and flight envelope are shown as a result of aerodynamic 
effects only; the effects of atmosphere on engine performance were not considered. The results show that the effects 
on aerodynamics are significant. Figure 3 shows trajectories for the four different atmospheric models analyzed: 
standard day, Edwards cold day, Edwards standard day, and Edwards hot day. The flight envelope was calculated 
for standard day and Edwards cold day only. 

The choice of atmospheric model had an effect on the optimal trajectory, especially in the high energy region, 
(fig. 3). At the higher energies, the maximum P s occurs at a higher altitude for the warmer days. The differences 
are due to the differences in densities at those altitudes between the atmospheric models. The temperature affects 
the density, which affects the calculated lift and drag. This is also reflected in the flight envelope, which shows the 
effects of atmospheric model on the dynamic pressure boundary, which is a function of the density of the air, as well 
as the thrust boundary at high energies. 

Figures 4(a), (b), and (c) show the effects of the different atmosphere models on the performance of the vehicle. 
In agreement with the data shown on figure 3, figure 4(a) shows that the vehicle must climb higher on warmer days 
in order to maximize its cost function. These higher altitude trajectories take more time and range (fig. 4(b)) and 
use more fuel (fig. 4(c)) than the lower altitude trajectories to achieve the same energy. They are, therefore, less 
efficient. If a large degree of accuracy is required in the calculations, the atmosphere must be modeled as close to 
the actual conditions as possible. 

Effects of Gravitational Models 

The effect of the two different gravitational models on flight envelope and trajectory is shown in figure 5. The 
two gravitational models are the altitude dependent model, based on a simple Newtonian force equation (eq. 34), 
and the simpler model in which gravity is assumed constant (cq. 35). Figure 5 shows that the flight envelope and 
optimum trajectory are overlaid for the two gravity models. 

Effects of Equations of Motion 

Figure 6 shows the effects of the differences in the equations of motion models on the flight envelope and optimal 
trajectory. For the most general equations, the angle between the thrust vector and the body axis was increased until 
a noticeable effect was seen on the trajectory calculations. This angle was determined to be approximately 6°. At 
angles up to 6°, there was only a small effect on the flight envelope and trajectory. The only difference is in the 
anglc-of-attack boundary of the flight envelope. The case with the thrust aligned with body axis overlays the case 
with the 6° angle between thrust and body axis, but both arc slightly different from the zero angle -of- attack case. 
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Figures 7(a) and (b) show that simplifying the equations of motion had little effect on the calculated performance 
of the vehicle. It took a thrust misalignment angle of 6° to show an effect on performance, in which case a slight 
increase in range and fuel used were needed to achieve the same energy state as the simpler eases. 

Effects of Constraint Conditions on Trajectory Calculations 


Constraint conditions had some effect on the trajectory, as figure 8(a) shows. For the constrained to level flight 
and the flightpath constrained eases, a condition of the constraints was that the optimization pick a point where the 
forces in the vertical direction arc balanced. Because of this, limiting the flightpath to realistic angles had the same 
effect as constraining the trajectory to level flight. The unconstrained ease, however, shows that a different optimal 
trajectory was picked. 

There is no significant difference in aircraft performance between the constrained to level flight case and the 
flightpath constrained case. However, both differed from the unconstrained case. The constraint conditions caused 
a less favorable but more realistic trajectory to be chosen. Figures 8(b) and (c) show that the constrained eases take 
more range and fuel to achieve flic same energy state than the unconstrained case does. 

CONCLUSIONS 

This paper examined the effects of various modeling assumptions on (light envelope generation and optimal 
trajectory estimation. The parameters investigated included atmospheric models, gravitational models, equations of 
motion, and constraint conditions. 

The different types of atmospheric models showed a definite effect on the trajectory, especially at high energies. 
Of the four examined, the colder days had better performance characteristics. The warmer day models took more 
time and distance to achieve the same energy state. If a large degree of accuracy is required in calculations, the 
atmosphere should be modeled as close to the actual conditions as possible. 

The two gravitational models examined were gravity as a function of altitude and a constant valued gravity force. 
They showed no difference in effect on flight envelope or trajectory calculations. 

Three different sets of equations of motion were examined. The assumption of zero angle of attack showed 
a noticeable but small effect on the anglc-of- attack boundary of the flight envelope, but showed no effect on the 
trajectory. For small angles between the thrust vector and the body axes, simplifying the equations of motion in the 
vertical plane had very little effect on calculated performance. 

The different constraint conditions had an effect on the trajectory. Constraining the trajectory to level flight 
limited the possible points to be chosen for the optimal path, and so made a difference in the optimal trajectory. The 
further constraint of limiting the flightpath to reasonable limits had no effect compared to the constrained to level 
flight ease. . — -- . 

Ames Research Center 

Dryden Flight Research Facility 

National Aeronautics and Space Administration 

Edwards, California, January 3, 1990 
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Figure 1. Sample flight envelope and optimal trajectory for a high-performance vehicle. 
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Figure 3. Effect of atmospheric model on flight envelope and trajectory. Flight envelope calculated for standard day 
and Edwards cold day only. 
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Figure 4. Effects of atmospheric model, 
(a) Altitude as a function of time. 
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Figure 4. Continued. 
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(c) Fuel used as a function of time. 
Figure 4. Concluded. 
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Figure 7. Concluded. 
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Figure 8. Effects of constraint conditions. 
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